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The r-modes of accreting neutron stars could be a detectable source of persistent gravitational 
waves if the bulk viscosity of the stellar matter can prevent a thermal runaway. This is possible 
if exotic particles such as hyperons are present in the core of the star. We compute bulk viscous 
damping rates and critical frequencies for r-modes of neutron stars containing hyperons in the 
framework of relativistic mean field theory. We combine the results of several previous calculations 
of the microphysics, include for the first time the effect of rotation, and explore the effects of 
various parameters on the viability of persistent gravitational wave emission. We find that persistent 
emission is quite robust, although it is disfavored in stars below 1.3-1.5 Mq depending on the 
equation of state. In some cases persistent emission is compatible with temperatures as low as 
10^ K, observed in some accreting neutron stars in quiescence. 

PACS numbers: 04.30.Db, 04.40.Dg, 26.60.+C, 97.10. Sj 


I. INTRODUCTION 

The r-modes (fluid oscillations governed by the Corio¬ 
lis force) of rapidly rotating neutron stars have attracted 
much interest as possible sources of gravitational waves 
and mechanisms for regulating the spins of neutron stars. 
See Ref. 0| for a recent review of the many physical and 
astrophysical issues related to the r-modes; here we focus 
on gravitational wave emission. Gravitational radiation 
drives the r-modes unstable and could lead to detectable 
gravitational wave emission in two scenarios. In one sce¬ 
nario, a newborn neutron star could radiate a substan¬ 
tial fraction of its rotational energy and angular momen¬ 
tum as gravitational waves changing in frequency on a 
timescale of a year or more Q. In the other scenario, 
r-modes in rapidly accreting neutron stars in low-mass 
x-ray binaries (LMXBs) could be persistent sources of 
periodic gravitational waves Hi (Here the r-modes 
provide a specific mechanism for a more general torque 
balance argument Hi) 

Currently the latter scenario for gravitational wave 
emission (from LMXBs) looks like a brighter prospect 
for detection. There is now evidence from several ap¬ 
proaches that the amplitude of an r-mode growing due 
to the instability is limited by nonlinear fluid dynam¬ 
ics to a relatively small value Ullillilllllllll 
While low-amplitude long-lived r-modes in newborn neu¬ 
tron stars still can lead to astrophysically interesting ef¬ 
fects such as the regulation of spins, a low mode ampli¬ 
tude renders the gravitational wave signal undetectable 
unless there is a very nearby supernova [3| • Also, if neu¬ 
tron stars contain particles more exotic than neutrons 
and protons—such as hyperons, where an up or down 
quark in a nucleon is replaced by a strange quark—there 
are additional viscous damping mechanisms which may 
eliminate the instability alto geth er in very young, hot 
neutron star s ll4L ll^ HR HTl HR H9l 1^ 1^ or strange 
quark stars [^. However, for the LMXB scenario, low 
amplitude is not a problem Q and the additional vis¬ 


cosity actually renders stars with exotic particles better 
candidates for gravitational wave detection. The reason 
for the latter is subtle, and bears explanation. 

Real neutron stars are not perfect fluids, and thus 
viscous (and other) damping mechanisms compete with 
gravitational wave driving of the r-modes. The strengths 
of the driving and damping mechanisms can be expressed 
as timescales which depend on the rotation frequency and 
temperature of the star (usually assumed to be nearly 
isothermal). Therefore it is useful to consider the location 
of a star in the temperature-frequency plane, as shown 
in Fig. n] The strength of viscosity can be graphically 
represented by a curve in that plane that is the locus of 
all points where the driving and damping timescales are 
equal—this defines a critical frequency as a function of 
temperature. Since the driving timescale decreases with 
frequency, stars above the critical frequency curve have 
one or more unstable r-modes, while stars below it are 
stable and stars on it are marginally stable. 

Figure n plots examples of critical frequency curves in 
the temperature range 10® K<C T ^ 10^° K appropri¬ 
ate for the LMXB emission scenario. (The needed range 
is higher than most observed temperatures of LMXBs 
in quiescence because of a thermal runaway; see below.) 
At low temperatures the damping in Fig.^is taken to be 
dominated by shear viscosity in a boundary layer between 
the solid crust and fluid core (from Ref. [^, augmented 
by a constant relative crust-core velocity in the range dis¬ 
cussed in Ref. |2^). There are many other possible curves 
for the low temperature part of the plot, corresponding 
to more complicated damping mechanisms such as turbu¬ 
lence or superfluid magnetoviscous effects [^, but 
they generally share the qualitative property of decreas¬ 
ing with temperature. (Some low-temperature curves lie 
above the observed range of spins entirely, but the ob¬ 
served spins of low-mass x-ray binaries are easier to ex¬ 
plain if the curve lies in the region indicated illlll.) At 
high temperatures the damping is probably dominated 
by bulk viscosity, either from the Urea process (pertur- 
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bation of /3-equilibrium) or from nonleptonic processes 
involving strange particles such as hyperons. The Urea 
process, which requires no exotic particles, does not af¬ 
fect the critical frequency curve for T < 10^° K and thus 
does not show up in Fig.^] Thus the top plot in Fig.Q] 
shows a critical frequency curve for a neutron star, and 
the bottom plot shows such a curve for a star with hyper¬ 
ons (with the high-temperature part of the curve derived 
from Ref. 0). This type of plot looks encouraging for 
gravitational wave detection because there is plenty of 
room above the curve for stars to be unstable and thus 
emitting gravitational waves. 

However, the gravitational wave emission duty cycle 
could be much smaller than 100% due to a thermal run¬ 
away 1^ . This happens generically when the criti¬ 
cal frequency decreases with temperature. In that case, 
plotted at the top of Fig. ^ a star will execute a loop as 
shown and radiate only during the time it spends above 
the curve. A stable star which begins at the bottom left 
of the loop is spun up by accretion until it moves above 
the critical frequency and the instability is triggered. The 
shear from the growing r-modes then causes the star to 
heat up, moving it rightwards on the loop. As the tem¬ 
perature rises, so does the rate of neutrino cooling, caus¬ 
ing the star to drop down toward the critical frequency 
again at a temperature of a few times 10® K. After falling 
below the critical frequency, the r-mode heating is re¬ 
moved and the star drifts leftward along the bottom of 
the loop until returning to its initial position. Although 
the time it takes for a star to complete the loop is dom¬ 
inated by the accretion rate, the timescale for gravita¬ 
tional wave emission depends mainly on the saturation 
amplitude of r-mode oscillations. For a saturation am¬ 
plitude (a in the notation of Ref. [^) of order unity, 
the duty cycle for gravitational wave emission is of the 
order 10“® [2^. Whereas the duty cycle can be as high 
as about 30% for the lowest predicted values of the satu¬ 
ration amplitude {a ~ 10“®), for typical estimates of the 
saturation amplitude {a ~ 10“^) the duty cycle is only of 
order 10“® |^. Advanced LIGO will be able to detect at 
most one LMXB (Sco X-1) without narrowbanding (and 
hurting its ability to see other sources), or 6-7 LMXBs 
with narrowbanding around a series of different frequen¬ 
cies 1^. The small number of detectable systems and 
the fact that the timescale for a star to complete a loop 
is much longer than a human lifetime mean that a duty 
cycle of order 10“^ or less is pessimistic for gravitational 
wave searches for the r-modes. 

If the critical frequency increases with temperature as 
in the bottom plot of Fig. ^ the thermal runaway can 
be blocked. A rapidly accreting star in an LMXB can 
have a duty cycle of order unity for emission of gravi¬ 
tational radiation as it sits on the curve or makes small 
peregrinations about it [s^. (It could also keep emitting, 
although less detectably, for some time after the rapid ac¬ 
cretion shuts off 13 .) This could happen for stars which 
exhibit a rise in the critical frequency curve at a low 
enough temperature. The rise is typical for stars where 
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FIG. 1: Critical frequency curves are given (qualitatively) 
by the dashed lines. The top plot includes no bulk viscosity 
due to hyperons or other strange particles. In this case an 
accreting neutron star traversing the loop indicated undergoes 
a thermal runaway and has a low gravitational radiation duty 
cycle. The bottom plot includes hyperon bulk viscosity. In 
this case the thermal runaway is blocked, and an accreting 
star is a source of persistent gravitational waves as it remains 
in equilibrium at the last arrowhead. 


high bulk viscosity processes involving hyperons 
or strange quarks are at work. Thus the case for the 
gravitational wave emission scenario is in fact strength¬ 
ened by high bulk viscosity from processes which are fun¬ 
damentally quark-quark interactions. One might expect 
similar critical frequency curves to arise for stars contain¬ 
ing other forms of strange matter such as a kaon conden¬ 
sate or mixed quark-baryon phase, but this has yet to 
be investigated. The thermal runaway is blocked if the 
increase in temperature (width of the loop in Fig. ^ is 
enough to take the star from the negatively sloped part of 
the instability curve to the positively sloped part. The 
larger the saturation amplitude of r-modes the greater 
is the increase in the temperature. Whether the star 
makes the jump to the positively sloped curve also de¬ 
pends on which cooling mechanisms are operative, and on 
the shape of the low-temperature instability curve, which 
in turn depends on which damping mechanism dominates 
at low temperatures. Notwithstanding the wide range of 
estimates for a, the shape of the negatively sloped curve, 
and the cooling mechanisms, the question of the r-modes 
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in LMXBs as a persistent source of gravitational radia¬ 
tion comes down to asking whether there is a rise in the 
instability curve around 10® K or lower. 

Our purpose in this paper is to revisit the question 
of whether there is a rise in the critical frequency curve 
around 10® K or lower, as in the bottom of Fig. ^ lead¬ 
ing to persistent gravitational wave emission. We focus 
(for now) on neutron stars containing hyperons (hyperon 
stars), because hyperons are in some sense the most con¬ 
servative and robust of the many proposals for exotic 
matter in the cores of neutron stars: Some properties 
of hyperons can be measured in the laboratory, both in 
vacuum and in the environment of a light nucleus, which 
combined with astronomical observations allows one to 
constrain some of the many uncertainties in building an 
equation of state Equations of state which allow 

for hyperons generally produce them at densities relevant 
for neutron stars, about twice nuclear density and up. 

We synthesize and extend results of previous work on 
this topic. In arguing the case for persistent gravita¬ 
tional wave emission. Wagoner and Reiseneg- 

ger and Bonacic base their critical frequency curves 
on bulk viscosity coefficients obtained by combining the 
results of Lindblom and Owen llll thereafter LO) and 
Haensel, Levenfish, and Yakovlev (hereafter HLY). 
The LO and HLY viscosities were obtained by differ¬ 
ent calculations and produced somewhat different re¬ 
sults. LO used a detailed self-consistent model of a multi- 
component fluid described by relativistic mean field the¬ 
ory , but since they were primarily concerned with the 
high-temperature regime appropriate to newborn neu¬ 
tron stars they treated the effect of superfluidity with a 
rough approximation. LO also made some errors which 
resulted in a bulk viscosity coefficient a factor of 20 or 
more too high. HLY were more careful with superfluidity, 
using consistent damping factors in the collision integrals, 
but instead of evaluating reaction rates and thermody¬ 
namic derivatives within relativistic mean field theory 
they used “order of magnitude estimates” of some quan¬ 
tities which resulted in more than an order of magnitude 
disagreement with the microscopic results of LO (and 
with Jones El). We correct some mistakes in LO and 
combine their self-consistent microphysical model with a 
more careful treatment of superfluidity similar to that of 
HLY. We also treat the macroscopic physics more care¬ 
fully than LO, including the effect of rotation on stel¬ 
lar structure which in some cases can significantly affect 
damping timescales and critical frequency curves. (HLY 
made only order of magnitude estimates of mode damp¬ 
ing timescales and did not plot critical frequency curves.) 

We also address the question “Does the viability of 
persistent gravitational wave emission require fine tun¬ 
ing of parameters?” Many of the parameters that go 
into building the equation of state have significant un¬ 
certainties, and those that go into computing reaction 
rates and bulk viscosities are even more uncertain. We 
investigate the viability of persistent gravitational wave 
emission with respect to variation of several microphysi¬ 


cal numbers such as hyperon coupling constants and the 
superfluid bandgap, and even which nonleptonic reaction 
is most important (ta king into account the results of van 
Dalen and Dieperink [IJ). We also investigate the de¬ 
pendence on the mass of the star, since cooling observa¬ 
tions E3 timing of radio pulsars in binaries E2l 
indicate a wider mass range than the traditionally as¬ 
sumed clustering around 1.4 Mq. (Cooling observations 
also might be interpreted to favor the existence of strange 
particles such as hyperons, although the data still can be 
fit by exotic cooling from purely nucleonic matter.) 

The organization of the rest of the paper is as fol¬ 
lows. In Sec. H we describe the microphysical model 
which leads to the equation of state and ultimately the 
macroscopic coefficient of bulk viscosity. In Sec. HI we 
plot the critical frequency curves for a range of neutron 
star masses and microphysical parameters. In Sec. IV 
we summarize our findings and discuss possible improve¬ 
ments. 


II. MICROPHYSICS 

We model the composition of a neutron star within 
the framework of relativistic mean field theory industri¬ 
alized by Glendenning and described in Ref. [33 ■ Here we 
do not consider meson condensates or exotica such as a 
mixed quark-baryon phase, but do allow for the presence 
of hyperons which are somewhat constrained by labora¬ 
tory and astronomical data. 


A. Equation of state 

There are several free parameters in this framework, 
some of which are better known than others. We use a 
range of parameters set by combining the old labora tory 
constraints with recent astronomical observations pI9|. 
At nuclear density (2-3x10^® g/cm®) the matter is, as in 
all models, composed of neutrons, protons, electrons, and 
muons. As the density rises hyperons generically appear, 
their order of appearance (and density thresholds) chang¬ 
ing somewhat with the precise set of parameter values 
but roughly corresponding to S“, A, E®, E+, 2“, start¬ 
ing at roughly twice nuclear density. Other baryons typ¬ 
ically appear in this framework only at densities higher 
than those found in the cores of the most massive neu¬ 
tron stars, nearly 10 times nuclear density. At such high 
densities the asymptotic freedom of the strong nuclear 
force is likely to lead to quark matter anyway, and any 
baryonic model must be considered somewhat suspect. 
As we shall see, most of the astrophysically important 
dissipation comes from lower densities (2-3 times nu¬ 
clear density) where this deficiency of the model is not 
important, i.e. low-energy effective field theory is good 
enough. The framework of relativistic mean field theory 
has some advantages over others, such as an inherently 
causal equation of state even at high densities, but has 
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some disadvantages including the neglect of correlations 
between particles which should become more important 
at high densities. 

In this framework, the strong interaction between 
baryons is modeled as a tree-level exchange of isoscalar 
mesons cr, vector mesons w, and isovector mesons p. The 
effective Lagrangian includes kinetic and tree-level inter¬ 
action terms for the baryons, leptons, and mesons, as well 
as an effective potential for the a expanded up to 
The variables of the theory are the Fermi momenta ki 
of the baryons and leptons, supplemented by the field 
strengths of the a, w, and p mesons. Expansion coeffi¬ 
cients and coupling constants can be fit to numbers ex¬ 
tracted from measurements of nuclei and hypernuclei at 
saturation density. The fitti ng p rocess (and Lagrangian) 
is described in detail in Ref. 133 . We use an updated set 
of fit parameters from Ref. . 

For N species of baryons, two of leptons, and three 
of mesons, the composition of the matter consists of 
TV -|- 5 unknowns (the baryon, lepton and meson fields) 
which must be determined by -I- 5 equations. The first 
three equations are the Euler-Lagrange equations for the 
mesons determined by varying the effective Lagrangian 
of Ref. 13 . The equations are rendered tractable by as¬ 
suming that the fields are given by their mean values in 
a uniform static ground state, and take the form 


Wo 


P03 

mla 



Ps,B km g fj c') 

B 


( 1 ) 

( 2 ) 

Cga{9aCrf-{^) 


Here the index B labels baryon species, wq is the timelike 
component of the vector meson field (the spatial compo¬ 
nents vanish), ns is the number density of baryon B, and 
P 03 is the isospin 3-component of the timelike component 
of the isovector meson field. Throughout this Section we 
use units such that h = c = 1. The scalar density is given 

by 


Ps,B = 


2Jb + 1 


r^B 




B 


27’"^ Jo \/k'^ + rn*g (a) 


k^dk. 


( 4 ) 


The coupling constants between the baryon B and the 
mesons are given by r , gujB, and gpB- We as¬ 
sume, as in Ref. |3li fbat they are given by one 
set of values {9(t,9uj,9p) for the nucleons and another 
{xcrgcr,Xuj 9 uj,Xpgp) for the hyperons, so that the x’s mea¬ 
sure relative coupling strengths for hyperons. The cou¬ 
pling constants for the self-interaction terms of the scalar 
field are b and c. The masses of the mesons are given by 
rria, and rup. The quantities for the Bth baryon 
I 3 B and Jb are, respectively, its 3-component of isospin, 
and its spin. The Fermi momentum of the baryon 
species is given by fcs while its effective mass is given 
by m*B = rriB — 9aBC!- The quantity rriN represents the 


average mass of a nucleon and is used to make b and c 
dimensionless. The constraints of charge neutrality and 
conservation of baryon number (if a particular baryon 
number density n is assumed) provide two more equa¬ 
tions which can be written, respectively, as 


'^mqi = 0 , 

i 

( 5 ) 

y^nB=n. 

( 6 ) 


B 


Here rii denotes the number density of fermion species 
i and qi is its electric charge in units of e. Now we are 
down to N unknowns, which can be eliminated with the 
N equations of generalized /3-equilibrium 

Pi = biPn - QiPe- (7) 

(The neutrinos are assumed to have zero chemical poten¬ 
tial.) Here pi is the chemical potential of fermion species 
i and bi is its baryon charge. (Actually there are -|- 2 
equations of /3-equilibrium including the leptons, but the 
equations for the electron and neutron are identities and 
thus do not count toward the elimination of unknowns.) 
The chemical potential of each baryon species is given by 


PB = 9u^BUJo + PpBPoshB + \Jk%+ m*^, ( 8 ) 

while that of each lepton species is given by 

PL = ijkl+ml- (9) 

Here k^ and rriL are, respectively, the Fermi momentum 
and mass of the Lth leptonic species. 

For a given baryon number density n, the baryon, me¬ 
son, and lepton fields are found by simultaneously solv¬ 
ing Eqs. 0 - 0 . The numerical technique is based on a 
multi-dimensional root finder and thus requires an initial 
guess to ensure finding the correct root. We first solve for 
a low value n below saturation density where the system 
can be approximated as a weakly interacting Fermi gas 
and analytic approximations can be found as in Ref. |38| . 
The code then steps up in baryon density n, at each step 
using the previous step’s values for the unknowns as the 
initial guess for the root finder. Once the fields are found 
for a given n, the mass-energy density e and pressure p 
can be determined: 

e = ^bmN {g^crf + {gcrCrf + 

1 rkB ! - 

— J k'^ + m*^k^dk 

+ / \jk'^ + m\k‘^dk, ( 10 ) 

P = -^bniN (Pauf - jcig^a)"^ - 
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TABLE I: Parameters for five relativistic mean field equations 
of state used in this paper. 


Name 

K (MeV) 

m*/mN 

Xa 

H3 

300 

0.70 

0.60 

H4 

300 

0.70 

0.72 

H5 

300 

0.80 

0.66 

H6 

240 

0.70 

0.67 

H7 

240 

0.80 

0.68 



( 11 ) 


Here the summations are over the baryon species (B) 
and the lepton species (L). These expressions can be 
combined to produce tabulations of p{e), the equation 
of state. 

There remains the problem of choosing the values of 
the constants {gafrua)^, {g^jraujY, b, c, x^, 

Xui, and Xp. Each choice of constants produces a dif¬ 
ferent equation of state with different hyperon thresh¬ 
old densities, neutron star maximum mass, etc. We use 
five equations of state from chosen to be com pat ible 
with recent measurements of neutron star masses 
and a gravitational redshift as well as hypernuclear 
data [ 33 ■ The constants in Eqs. 0-0 are fit in the 
manner of Glendenning to a A binding of -28 MeV, 
saturation density 0.153 fm“^, binding energy per nu¬ 
cleon -16.3 MeV at saturation, and isospin asymmetry 
coefficient 32.5 MeV, as well as a range of incompressibil¬ 
ities K, nucleon effective masses m*, and hyperon cou¬ 
plings Xa-- Since it is easiest to qualitatively understand 
the equations of state in terms of the latter three parame¬ 
ters, we give them in TableHl Low K and high m* lead to 
a soft equation of state, while high K and low m* lead to 
a stiff equation of state. Since most of the matter in the 
star is at high density, the hard core repulsion dominates 
the strong interaction. Thus high x^ postpones hyperon 
formation to higher densities, reducing hyperon popula¬ 
tions (and thus viscosity) and stiffening the equation of 
state. 

All of these equations of state are somewhat different 
from that used by LO and HLY. Though both HLY and 
LO use the “case 2” equation of state from an older paper 
by Glendenning |^. there were misprints in that paper 
which were corrected by HLY but not by LO. The scalar 
self interaction coupling constants b and c in Ref. 
should be 3 and 4 times larger, respectively, than the 
values that are quoted for them. Also, case 2 of Ref. M 
used K = 285 MeV and x^j = x^j = Xp = y^2/3, and 
had several minor differences in other constants (we use 


values from the 2002 Particle Data Group review 
As a result, the newer equations of state we use have 
higher hyperon populations and thus higher viscosities. 


B. Relaxation timescale 

The reactions that contribute most to bulk viscosity 
are the weak interaction processes, as their relaxation 
timescales are within a few orders of magnitude of the 
r-mode period (milliseconds). 

Among these, the most significant that can be calcu¬ 
lated from first principles are the non-leptonic weak in¬ 
teractions involving the lightest hyperons, and A, as 
they occur at lower densities than the more massive par¬ 
ticles and thus have higher populations in a given star. 
Following LO, we calculate matrix elements for the reac¬ 
tions 


n + n p -I- S (12) 

n+p ^ p + K (13) 

as tree-level Feynman diagrams involving the exchange of 
a W boson. The latter reaction was also considered (with 
a phenomenological free parameters) by HLY. These re¬ 
actions are combined to calculate an overall microscopic 
relaxation timescale. (This is the timescale on which a 
small perturbation of the neutron fraction returns to its 
equilibrium value.) 

In addition to these reactions, Jones includes 

n -|- n n + A (14) 

which is the dominant nonleptonic process observed in A 
hypernuclei in the laboratory. We do not consider this 
reaction as it has no contribution based on a W-boson 
exchange, even though the bare-mass interaction rate for 
this process is known. Using Eq. (15) from Jones 
we estimate that reaction da has a relaxation timescale 
that is longer than the timescale of the A reaction m 
by a factor of about 2 and thus does not change the 
final viscosity much. (When accounting for the macro¬ 
scopic structure of the star, the S” process dominates the 
overall bulk viscosity since the population extends to 
lower densities and therefore a greater volume fraction of 
the star than the A population.) There are various other 
concurrently occuring processes which contribute to the 
net reaction rate but their rates are not easy to predict. 
Even these rates are subject to substantial uncertain¬ 
ties: Van Dalen and Dieperink use a meson-exchange 
model for all three hyperon reactions and find reaction 
rates of order 10-100 times greater than LO and thus 
bulk viscosity coefficients 10-100 times lower than LO. 
(As we shall see in the next section, such a change has 
less of an effect on the r-mode critical frequency than one 
might think.) Our result for the net rate is then a lower 
limit, and consequently an upper limit on the bulk vis¬ 
cosity, in the typical range of temperatures for LMXBs. 
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To summarize the results of LO, the relaxation 
timescale can be computed as a function of the equi¬ 
librium matter fields. Neglecting superfluidity, the result 
is 


1 

T 


1927r3 




5^ 

nxn 


(15) 


when the A hyperons have not yet appeared and 


when both A and are present. Here kT is Boltz¬ 
mann’s constant times the temperature. The matrix ele¬ 
ments M. are computed as tree-level Feynman diagrams, 
squared, summed over initial spinors, and averaged over 
the angular part of the collision integral (indicated by 
0). (Note that van Dalen and Dieperink do not an¬ 
gle average.) They are functions of the Fermi momenta 
and meson fields given by Eqs. (4.28) and (4.29) of LO, 


1 

T 


jkTf 

1927r3 


(A:e(|A4||) -I- A:a(|A4a|)) —^ 


(16) 


J 


(|A4a|^) = 1120 (l - gl^) (l - ml_rnfm\ - 20 (l - gl^) (l -f gif,) mlm*p (S/Xp/rA - kl) 

-10 (1 + glp) (1 - glfl) m*pml (6/r„/ip - 3kl + kl)+2 [(l + gl^) (l + gif,) + 4g„pgpA] 

X [5/ip/iA — kf,) + kf, (lO/ira/Xp -|- 5kl + lO/Cp — ^a)] + [(l + 5np) (l + 5 pa) ~ ^gnpOpA] 

X [lOgntJ-A -b /ci) -b k\ (-20/^^ -b 15fcp - ik\ + b{kl - klf jikl - kj,))] } , (17) 

(|A4e|^) = {l80 (l — glp) (l — gl-^) — 40 (l — glp) (l -b gl^) (Slin/rs — A:|) — 20 

^ (l + glp) (l ~ (6^nMp — 3fcp -b /c|) — 5 (l — g^p) (l — ^p™E — 3fcp — A:|) 

+4 [(l + glp) (l + ffns) + 4gnp5nE] [lO/T^ + fcl) + (6fc^ — Ski ~ (lO^n + ^^1 

-ks)] + [(1 + glp) (1 + gls) - 45'„p5nE] [-lOmf {Sgp^Y + k%) -b 10ir„ {Q^ngp^Y - ‘^gpk^ - SgYkl 

+Ms^x:) + IS/Cpfcl — Sk^) } . (18) 


We reproduce them here since the latter equation in LO is 
missing some terms (although the code used to generate 
the LO results is not missing them). Also, LO actually 
used kinetic energies instead of chemical potentials which 
should be used to be consistent with the quasiparticle 
picture of Fermi liquid theory. Correcting this reduces 
the relaxation time by a factor of 10 or more, making 
up much of the difference with the results of van Dalen 
and Dieperink who use a drastically different model 
of the interaction. In the matrix elements, Gp is the 
Fermi constant and 9c is the Cabibbo angle. The g’s are 
axial-vector coupling constants, which we take to have 
their values measured from vacuum /3-decay of particles 
at rest ^3. The remaining factor on the far right of 
Eqs. 03 and HU) is explained in the next subsection. 

HLY, in contrast, used the nonrelativistic limit 

(|MsP) =4G|sin2 20c(l + 3g„p5„E)" (19) 

of Eq. (I18II and left the term in parentheses as a phe¬ 
nomenological parameter, which they simply set to 0 . 1 . 
This is somewhat justified because the Particle Data 
group values for the couplings gnp and result in a 
near-perfect cancellation of the term. Any in-medium 
change of these values would have a disproportionately 
large effect as a result. We (and LO) find that including 
the full Eq. (I18II erases the effect of the near-cancellation 


of the leading-order terms and produces a running value 
of the term which can increase by more than an order of 
magnitude above 0 . 1 . 


C. Bulk viscosity 

We use for C, the macroscopic coefficient of bulk vis¬ 
cosity, the special relativistic expression derived by Lind- 
blom and Owen [l^ : 


> ^ P{loG - 1 o)t 
1 -b (wr)^ 


( 20 ) 


where p is the pressure, uj is the angular frequency of the 
r-mode in a frame corotating with the star, and r is the 
net microscopic relaxation time we just computed. The 
fast adiabatic index 700 (for infinite-frequency perturba¬ 
tions) can be written as 



( 21 ) 


where rii = kf/{3Tr'^) is the number density of species 
i and the partial derivatives can be evaluated explicitly 
from Eq. HU)- The slow adiabatic index 70 (for zero- 
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frequency perturbations) can be written as 


D. Superfluidity 



and is straightforward to evaluate for example by dif¬ 
ferentiating the equation of state. However, while LO 
derived the correct relativistic expression and noted that 
it reduces the factor — 7 o by a factor of 2-2.5 from 
the nonrelativistic expression, their Fig. 2 and code for 
the rest of the paper incorrectly used a data file with the 
nonrelativistic expression and thus are a factor 2-2.5 too 
high. 

Now we address the factor Sfi/Sxn in Eqs. llT^ and itTCI) 
in the same manner as LO. This factor is determined 
by the constraints of charge and baryon conservation 
and ®, plus the constraint that the reaction 

n -b A p"*" -b (23) 

is in equilibrium since its reaction rate is many orders 
of magnitude greater than the weak interaction rates. 
This implies that both the non-leptonic reactions have 
the same chemical potential imbalance: 


6 fi = 5nn - SfiA = 26fj,n - SfXp - Sfi^. (24) 

Assuming small perturbations and using the neutron 
fraction this yields the relation 


6fj, 


nsSxr, 


- Qlri 


^An 


(/^n PA){_C)^np C^Ap 

2(3A - /?P - /?E 

Pp /?I])(Q;nA O^Aa) 
2 /?a — Pp — Pi] 


, (25) 


where f3i is given by 

— OLjii “b OiAi (26) 

and Ui is given by 


OLij 



(27) 


These expressions apply for densities where the A and S“ 
hyperons are both present. In the regime where only the 
hyperon is present the strong interaction constraint 
is no longer a factor and one obtains a simpler form 


2(5/i 

TlAiSXn 


dcTnn 2(^CXpfi “b O^ln “b CXjip F 

Fopp + asp + OpE + asE. ( 28 ) 


LO used numerical differencing to calculate the deriva¬ 
tives aij, but due to a coding error the a meson field 
was not differenced properly. The maximum error this 
induces in the value of 6 fj,lnB6xn is less than 0.3%, well 
below the uncertainties of the problem. 


In the temperature range of interest (below 10^° K), 
nucleons and hyperons are expected to form Cooper pairs 
near their Fermi surfaces and act as superfluids. This 
greatly slows reaction rates and has an important effect 
on transport coefficients such as bulk viscosity. 

The energy associated with the Cooper pairing is 
given by the bandgap. Assuming ^Sq pairing, the zero- 
temperature bandgap A n is related to the superfluid crit¬ 
ical temperature Tq by jig 

kTc = 0.57Ao. (29) 


We use the LO fit to the finite-temperature bandgap 


A(T) = Ao 



(30) 


For the A hyperons we use the empirical fit made by 
LO to the zero-temperature gap function Aqa of the A 
hyperon as computed by Balberg and Barnea ^3- The 
calculation by Balberg and Barnea is constrained by 
experiments on double A hypernuclei, though like other 
calculations of Aa it is likely only good to within a factor 
of 2 or 3. The zero-temperature gap depends on the total 
baryon number density n and on the Fermi momentum 
kA in a way that LO found was well fit by 


AoA{kA,n) = S.lfcA^ (1.52 -/ca)^ 


0.77 F 0.043 (6.2n-0.88)^ .(31) 


Here Aqa is in MeV, kA in fm“^, and n in fm“^. Thus the 
critical temperature for A superfluidity peaks somewhat 
below 10^° K. 

The superfluid bandgaps are not as well known as 
Aa due to the absence of similar experiments on hyper¬ 
nuclei containing the hyperon. Determination of the 
hypernuclei energy levels for the S“ hyperon are con¬ 
strained by the very short decay time of the hyperon. 
Takatsuka et al ^3 have however calculated, using sev¬ 
eral models of the nuclear interaction, that the bandgap 
lies in the range Aa < A^- < IOAa. To account for 
the uncertainty in the bandgap we perform our calcula¬ 
tion for two cases as in LO: when Ae;- = Aa and when 
Ae- = IOAa. 

For the superfluidity of protons we use model 2p of 
Ref. which is one of the bandgap models used by 
HLY. (LO did not model proton superfluidity, since it 
was not relevant at the temperatures greater than 10^° K 
that they considered.) This bandgap model is not based 
on a specific calculation of the bandgap but preserves the 
general features of bandgaps predicted by various micro¬ 
scopic theories. In this model the zero temperature gap 
depends only on the Fermi momentum kp of the proton. 
The relation for Tc is given by 


Tc = To 


kl 

{kp - k2Y 

F kl 

_{kp - k2Y F 


(32) 
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where fci = 1.117 fm“^, k 2 = 1.329 fm“^, fca = 
0.1179 fm“^ and Tq = 17 x lO^K. The Fermi momentum 
kp is in fm“^ units and Tc is non-zero only if kp < ki. 
This leads to Tc for the protons peaking somewhat below 
10^° K, similar to the A critical temperature. 

Since neutrons are, after all, more abundant than other 
particles in neutron stars, they have higher energies and 
at supernuclear densities are unlikely to form ^5'o pairs. 
They can have a ^P 2 pairing as considered in Ref. [i^ . 
but it is much weaker and corresponds to a maximum 
critical temperature Tc ~ 3 x 10® K. Below this tempera¬ 
ture the contribution to the critical frequency curve from 
hyperon bulk viscosity is quite low, most of the contri¬ 
bution coming from viscosity at the crust-core interface. 
Therefore neutron superfluidity does not much affect the 
viability of persistent gravitational wave emission from 
the r-modes, and we neglect it here. 

The effect of superfluidity on the microscopic relax¬ 
ation timescale or (Cni appears, after doing the 

lengthy collision integrals (see HLY and LO), as mul¬ 
tiplicative factors in the average matrix elements: 

(lAdsP) ^ Rnnpi:{\Ms\'^), (33) 

(IMaH ^ RnppA{\MA\^). (34) 

If none of the participating particle species is superfluid, 
an R factor is 1. Each R factor is reduced as more 
species become superfluid. Our treatment of these fac¬ 
tors is based on that used by HLY, who calculate R for 
the cases when either one or two of the particles partic¬ 
ipating in the reaction are superfluid. HLY do not how¬ 
ever provide a 3-particle R which is required for calculat¬ 
ing the timescale of the A reaction m- Our numerical 
investigations found (Him that expressing RnnpY. as a 
product of two single-particle factors RpR^- instead of 
HLY’s two-particle factor does not affect the r-mode crit¬ 
ical frequency much. Therefore throughout this paper 
we express all i?’s as products of the appropriate single¬ 
particle factors. The expression for the single-particle R 
from Eq. (30) of HLY is 

R = -^-exp (^0.5068 - v^O.50682 -h , (35) 

where y = A{T)/kT, a = 1 -I- 0.3118y^ and 5=1-1- 
2.566y^. LO used a simple factor e~^ for the hyperons 
only, since they were focusing on the ordinary fluid case. 
The exponential is indeed the dominant behavior of the 
HLY expression for y ^ 1, but the HLY expression is sig¬ 
nificantly more accurate for y less than about 5-10, which 
turns out to be the astrophysically interesting region. 

We plot the most important results of our micro¬ 
physics calculations. Figure [3 shows a typical relax¬ 
ation timescale as a function of density, and Fig.|3shows 
the corresponding bulk viscosity coefficient. Among the 
equations of state we use, H3 has a relaxation timescale in 
the middle of the distribution. Neglecting superfluidity, 
the timescales are about a factor of 40 shorter than in LO. 
About 20 of this comes from making the matrix elements 



FIG. 2: Microscopic relaxation timescale as a function of 
density for various temperatures. The equation of state is H3 
and the and A bandgaps are assumed equal. Thus at the 
highest temperature plotted, no superfluid effects are present. 



FIG. 3: Bulk viscosity coefficient as a function of density for 
various temperatures, for the same model as Fig. 13 The star 
is rotating with angular frequency 2300 radians per second. 


consistent with the quasiparticle picture, and most of the 
rest comes from changes to the equation of state (includ¬ 
ing accounting for the typos in Ref. M which propagated 
to LO). Fig. 13 is reduced from LO by a further factor of 
2 due to the relativistic factor joa — 7o • Including super¬ 
fluidity, the shapes of the curves are different from LO 
but the peak timescales are comparable. It is difficult 
to compare to HLY because at these temperatures the 
high-frequency approximation they use (wr 1) does 
not hold. 


III. MACROPHYSICS 

With the bulk viscosity coefficient ( in hand, we can 
proceed to compute r-mode damping times and critical 
frequencies. For this we also need the hydrodynamic 
structure of the r-mode as well as equilibrium models of 
the structure of a neutron star (i.e. its density profile). 
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A. Stellar structure 


We are motivated by the sensitivity of the hyperon 
population to the central density to include the effects 
of general relativity and rotation in our stellar models. 
Relativistic stellar models have higher densities than non- 
relativistic ones, and rotating models have lower densi¬ 
ties than nonrotating ones. LO used relativistic, non¬ 
rotating stellar models; while other work such as HLY 
and Ref. did not treat the stellar or mode structure 
in much detail. Since we are interested in exploring a 
broader parameter space, we do not restrict ourselves to 
nonrotating stars. 

We take into account the effect of rotation on the neu¬ 
tron star structure using Hartle’s slow-rotation approx¬ 
imation 1^. Hartle’s formalism is based on treating a 
rotating star as a perturbation on a non-rotating star. 
We start by solving the Oppenheimer-Volkoff equations 
for a static spherically symmetric star in a form due to 
Lindblom 


dm ATTe{h)r{h)^[r{h) — 2m{h)] 

dh m{h) + ATrr{h)^p{h) 

dr r{h)[r{h) — 2m{h)] 

dh m{h) + 4TTr(h)^p{h)’ 


(36) 

(37) 


using as independent variable the specific enthalpy 


h{p)=f dp'/[p'+ e{p')]. (38) 

Jo 


From r{h), the radius at which the enthalpy is h, and 
m{h), the mass contained within a sphere of that radius, 
we find m{r). We also find ^(r), the logarithm of the 
time-time metric component which is based on h(r). We 
then use this solution of the static star to solve for a 
rotating star with the same central pressure. 

The metric of a slowly rotating star to second order in 
n, the spin frequency, can be written as 

ds^ = - e" [l + 2{hQ + h2P2)\dt^ 


2(mg-bm2P2)’ 

2M(r)' 

(r — 2M(r)) 

r 


+ [1 -|- 2(v 2 — h2)P2\ [dO^ + sin^ 6{d(j) — wdt)^] 

-b 0{n^). (39) 


and the boundary conditions that at the center of the star 
a;(0) = (Dj, an arbitrary constant, and {du){r)/dr)r=o = 0. 
From Eq. (P|l one determines fl and the angular momen¬ 
tum J corresponding to Cui as follows: 


6 \ dr 


J 

n = lu{R) + 


duj(r) 


r—R 


2J 


(42) 

(43) 


Here R is the radius of the star. To obtain the desired 
value of n one can scale uj thus: 


new /^old)- (44) 


The angular momentum J corresponding to ^new can 
then be calculated using Eq. 621). 

A measure of the change in the star structure caused 
by rotation is ^(r, 9) which represents the displacement 
of constant energy density surfaces between the rotating 
star and the corresponding non-rotating star. This dis¬ 
placement’s angular dependence can be separated out: 

?(A0)=Co(r)+6(r)P2. (45) 


The I = 0 part of this deformation, can be obtained 
by solving the following equations: 


dmo{r) 

dr 


dPojr) 

dr 


+ 

-b 


-b 

-b 


47rr2^|^ Hr) + p{r)]pl{r) 
1 -2/ ^ i (d}{r) 

1 odH{r) 

? —“ <’■>' 
mg [l -b 87rr^p(r)] 
r — 2M (r) 

47r [e(r) -b p(r)\ r^ 
r — 2M(r) 


Pair) 


r'^Hir) / dw{r)\ 
12[r —2M(r)] \ dr ) 
Id/ r^H{r)Y{r) \ 

3 dr \ r — 2 M{r) J ’ 


(46) 


(47) 


Here P 2 = P2(cos9) is the 2nd order Lagrangian poly¬ 
nomial. The quantity uj is the frame dragging frequency 
that is proportional to H, while hg, / 12 , mg, m 2 , and 1^2 
are all functions of r that are proportional to . Before 
we solve the Einstein equations for all the subscripted 
quantities in Eq. 63) we need to solve for w = fl — uj. 
This is done by solving 


1 d 
r^ dr 



du!{r) \ 
dr ) 


4 

H— 
r 


djjr) 

dr 


uj{r) = 0, 


(40) 


with the definition 


j(r) = exp 



2M(r) 


1/2 


r 


(41) 


where pg is a pressure perturbation, the density profile 
e(r) is the energy density of the non-rotating star as func¬ 
tion of r, and p(r) is the pressure of the star as a function 
of r. The boundary conditions for the above equations 
are that both mg and pg vanish at the origin. 

The gravitational mass of the rotating star is given by 

72 

M'{R) = M{R) + mg{R) + (48) 

which is greater than the mass of the nonrotating star 
with the same central pressure. We wish to build se¬ 
quences of stars rotating at various rates all with the 
same gravitational mass. To do this we first solve the 
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Oppenheimer-Volkoff equations; then, solve Eq. (EHll . 

and lITTI) : then calculate the new mass of the ro¬ 
tating star using Eq. (1^ . This procedure is repeated 
with lower values of the central pressure until the same 
mass as the nonrotating star is obtained. 

To perform the viscosity calculations we need informa¬ 
tion on the structure of the rotating star contained in its 
density profile. Once the correct mass is obtained the 
density profile e(r) is calculated by using 

^rotir) = €,tat{r) - (49) 

dr 

where 


?o(0 


-Pair) 


e{r) +p{r) 
dp{r)/dr 


(50) 


Here, erotif) is the density profile of the rotating star, 
while Cstatir) is the density profile of the non-rotating 
one. This yields only the spherically symmetric part of 
the density perturbation, which is the dominant one for 
our purposes. 


B. Driving and damping timescales 


The stability of an r-mode is determined by calculat¬ 
ing the damping and driving timescales. The mode is 
unstable if the driving timescale is shorter than the vis¬ 
cous damping timescale. This can be expressed in terms 
of an overall r-mode timescale Tr such that 
1 1 1 
Tri^l,T)~ TGRin)^ Tv{n,T)^ ^ ’ 

where tgr < 0 is the damping timescale of the mode 
due to gravitational radiation and > 0 is the damp¬ 
ing timescale resulting from the sum of all viscous pro¬ 
cesses in the star. The star’s spin angular frequency and 
temperature are represented by H and T, respectively. 
A mode unstable to gravitational radiation corresponds 
to Tj. < 0 in Eq. and the critical frequency flc of 
the star as a function of temperature is found by solving 
l/Tr{nc,T) = 0. 

The gravitational radiation timescale is given by 


1 _ if dE\ 


(52) 


where E is the mode energy in a frame rotating with the 
star and {dE/ dt)GR is the rate at which energy is emitted 
in gravitational waves in the same frame. For the cur¬ 
rent quadrupole r-mode, there exist simple expressions 
(to lowest order in mode amplitude and angular veloc¬ 
ity of the s tar) for the energy and gravitational radiation 
timescale: 


E 


1 

Tgr 


2 

1310727r 


/■« g 

/ e{r)r dr, 

Jo 

e(r)r® dr. 



(53) 


Here a is a dimensionless amplitude coefficient which 
cancels out of the driving and damping timescales. 

The viscous damping time scale Ty is given by 


1 _ 1 
Tv 2E 



(54) 


where {dE/dt)v is the rate at which energy is being 
drained from the mode by viscosity. Here we consider 
only hyperon bulk viscosity from the nonleptonic pro¬ 
cesses discussed in the previous section. Leptonic (Urea) 
processes also contribute, but in LO they were shown to 
be overwhelmed by the viscosity due to nonleptonic pro¬ 
cesses. Although the crust-core viscosity is important at 
low temperatures—it determines the shape of the nega¬ 
tively sloped part of the instability curve which in turn 
puts a constraint on how much farther along the tempera¬ 
ture axis the rising part of the instability curve should be 
located at to prevent a thermal runaway—we neglect it. 
For a reasonable estimate of the shape of the crust-core 
viscosity curve and of the value of the saturation am¬ 
plitude the stars temperatures hould increase to around 
10® K 1^. To explore the effect of the crust -core vis¬ 
cosity and the saturation amplitude in more details is 
beyond the scope of this paper. We can then compute 
the viscous damping timescale as in LO: To lowest order 
in O, we can write 


=—Att / ^(e(r))(|V • ^Ep)r® dr, (55) 

Jo 

where (| V ■ dw|®) is the angle average of the square of the 
hydrodynamic expansion. LO found that the expansion 
found numerically in Ref. could be fit well by 



(|V • 6 v\^) 



(56) 

where e is the mean density of the (nonrotating) star. 
[In Eq. (6.6) of LO the last factor is missing due to a 
misprint, but is properly included in the calculations.] 
This approximation breaks down for r/R ^ 1, but since 
the bulk viscosity is coming from the hyperon core that is 
enough. In Table HD we give typical radii of nonrotating 
stars and their hyperon cores. 


C. Critical frequency curves 

We are now in a position to solve for the critical fre¬ 
quency as a function of temperature for various neutron 
star models. This must be done numerically. 

Because of the inclusion of rotation the process is some¬ 
what more involved than in LO. The critical frequency 
curve is calculated in two stages to ensure that the correct 
gravitational mass is retained after including the correc¬ 
tions to it from rotation. In the first stage we calculate 


164 025 
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EOS 

Ec (10^'‘g/cm®) 

R (km) 

Rs (km) 

Ra (km) 

H3 

6.86 

13.74 

6.38 

0.00 

H4 

6.12 

13.85 

4.50 

0.00 

H5 

7.70 

13.33 

5.37 

0.00 

H6 

7.84 

13.44 

6.70 

2.46 

H7 

9.00 

13.01 

6.03 

0.00 


TABLE II: Parameters of non-rotating IAMq equilibrium 
stars for each equation of state (EOS). The central mass den¬ 
sity is Ec, R is the stellar radius, Rt. is the distance out to 
which E“ hyperons appear, and i?A is the corresponding dis¬ 
tance for A hyperons. 



FIG. 4: Critical freqnency curves I.IMq neutron stars with 
equal bandgaps for and A hyperons. The dashed curve is 
for equation of state H3 and the solid curve is for H6. The 
thinner curves neglect the effect of rotation on stellar struc¬ 
ture while the thicker curves include it. The curves including 
rotation are extended to their peaks, but no further. 


the critical frequency as a function of temperature for a 
nonrotating star using the Oppenheimer-Volkoff formal¬ 
ism. Then a set of points on this curve is taken. In the 
second stage, for each point we find the structure (density 
profile) of a star rotating at that frequency with the same 
gravitational mass as the nonrotating star. Since Har- 
tle’s formalism produces stars with greater gravitational 
masses, this step requires several iterative trials of the 
central pressure to find the star with the same gravita¬ 
tional mass. Using the new density profile with rotation, 
we keep the frequency fixed and use the corresponding 
temperature on the nonrotating curve as an initial guess 
for a routine that finds the root of 1 /r,., thereby obtaining 
the correct temperature for the marginally star rotating 
at that frequency. Because of this it is cumbersome to 
trace out parts of the curves where the critical frequency 
decreases with temperature, and since those parts of the 
bulk viscosity curves are not important for determining 
the viability of persistent gravitational wave emission we 
neglect them. 

The importance of rotation is shown in the critical fre¬ 
quency curves plotted in Fig. (Kepler frequencies for 
these equations of state are 750-850 Hz for 1.4 Mq stars.) 


Equation of state H6 is soft and has the highest hyperon 
population for a given mass, and thus its curve is higher 
than for H3 which is intermediate. This plot shows the 
general trend that the effect of rotation is more significant 
for softer equations of state. For persistent gravitational 
wave emission the key question is whether or not a neu¬ 
tron star undergoing thermal runaway will be blocked 
by hitting the curve at a temperature lower than a few 
times 10® K. Even for the soft H6 equation of state, the 
effect of rotation on stellar structure changes the temper¬ 
ature at which a star hits the curve by no more than a 
factor of 2 (still within the region that is good for gravi¬ 
tational wave emission). More significant is the fact that 
the frequency at which the curve peaks can be consider¬ 
ably reduced from the value neglecting rotation. If the 
peak is low enough, thermal runaway can still occur— 
for example in this plot it can occur for H3 stars if they 
hit above 700 Hz—and the outlook is bad for persistent 
gravitational wave emission. 

In Fig. 0 we show the effect of the equation of state 
on the critical frequency curve for 1.4 Mq stars. The 
curves generically reach their peak at temperatures of 1- 
2x10® K, which is good for persistent gravitational wave 
emission. However, the heights (frequencies) of the peaks 
change significantly. The peak is highest for H6, which 
has the largest hyperon population. It is almost iden¬ 
tical for H3 and H7, although those two equations of 
state are the furthest separated in the {K,m*,Xa) pa¬ 
rameter space. This seems to be due to the similar radii 
of their hyperon-containing cores as in Tabled All three 
of these equations of state then allow persistent gravita¬ 
tional wave emission in 1.4 Mq stars. The peak is lowest 
for H4, which is not surprising since that is the stiffest 
equation of state and has the lowest hyperon population. 
For H4 the peak is so low that, depending on the details 
of the crust-related damping, the full curve including that 
damping might show no significant effect from hyperon 
bulk viscosity—it might decrease throughout the tem¬ 
perature range and never allow persistent emission. In 
any case, for H4 and H5 persistent emission is not pos¬ 
sible for the most rapidly rotating stars in LMXBs (up 
to 619 Hz) if they are 1.4 Mq. This demonstrates in 
principle how observation of r-mode gravitational waves 
from an LMXB could rule out some equations of state 
or, more broadly, constrain hyperon populations. 

In Fig.Elwe plot critical frequency curves for stars with 
different masses for the H4 equation of state, which is the 
stiffest of those we use. (For other equations of state, the 
stellar mass affects the curve more but is constrained to 
a narrower range, so the overall variation is greatest for 
H4.) We also vary the bandgap from the A value 
to 10 times that value. It can be seen that increasing 
the mass pushes the curve to lower temperatures and 
raises its peak. This is due to the fact that higher mass 
stars have higher central densities and hence, in abso¬ 
lute terms, larger hyperon populations than their lower 
mass counterparts. Persistent r-mode gravitational wave 
emission is viable in all but the 1.4 Mq neutron stars. 
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FIG. 5: Critical frequency as a function of temperature for 
1.4 Mq stars with different equations of state. All curves 
assume the superfluid bandgaps of the A and E“ hyperons 
are equal. 



Temperature (K) 


FIG. 6 : Instability curves for H4 stars with different masses. 
The thicker curves use the larger E“ bandgap while the thin¬ 
ner curves use the smaller one. The solid lines represent 
2.0Mq stars, the dot-dash lines repsents 1.8Mq stars, the 
dotted line represents I. 6 M 0 stars, and the dashed line rep¬ 
resents I. 4 M 0 stars. 


The temperature at which this can happen is strongly 
dependent on the mass, and to a lesser but still signifi¬ 
cant extent (factor 3) on the E“ bandgap, with higher 
bandgaps leading to higher temperatures at a given fre¬ 
quency. (The small horizontal sections of the IAMq and 
I. 6 M 0 large bandgap curves should actually dip slightly. 
This behavior is missed due to our method of including 
rotation, but is unimportant for the thermal runaway and 
persistent gravitational wave emission.) 

The curves for 1.8-2.0 Mq stars are particularly in¬ 
teresting in light of predicted interior temperatures of 
neutron stars in LMXBs. These can be determined from 
observations in quiescence of the surface temperatures of 
three slowly accreting neutron stars mentioned in Brown, 
Bildsten, and Chang |^. To calculate the interior tem¬ 
peratures we use Fig. 2 from Yakovlev et al |5^, which 
depicts the relationship between the core and the mea¬ 
sured surface temperatures using different atmospheric 


FIG. 7: Variations on the critical frequency curve for a 
1.4 Mq H3 star with equal A and E~ bandgaps. The solid 
curve is our standard model. The dotted curve uses asymp¬ 
totic values for the axial couplings instead of vacuum at-rest 
values. The dot-dash curve uses the 2-particle superfluid re¬ 
duction factor for the E“ interaction instead of two 1-particle 
factors. The dashed curve simulates the fast timescales pre¬ 
dicted by van Dalen and Dieperink. 


and equation of state models. We thereby obtain a tem¬ 
perature range (1.5 x 10^-4 x 10®K) which covers LMXBs 
whose thermal emission spectra have been observed. This 
result is probably statistically biased, since LMXBs with 
higher accretion rates are likely to have higher tempera¬ 
tures but surface temperatures are harder to observe be¬ 
cause of the accretion discs. At any rate, the most mas¬ 
sive H4 stars have critical frequency curves which would 
allow gravitational wave emission scenarios 0113 to op¬ 
erate in that temperature range. Less massive stars allow 
it at higher temperatures, which may be consistent with 
the more rapidly accreting neutron stars such as Sco X -1 
whose temperatures are poorly known. 

The effects of varying some additional microphysical 
parameters are shown in Fig.[7| which depicts curves for 
1.4 Mq neutron stars with the H3 equation of state and 
the superfluidity bandgap equal to that of the A. 
The variations include using the asymptotic (quark) val¬ 
ues for the axial-vector couplings in the reaction rates, 
using HLY’s two-particle superfluid reduction factor in¬ 
stead of the product of two one-particle factors, and us¬ 
ing relaxation times an order of magnitude faster from 
van Dalen and Dieperink [IJ. All of these curves peak 
at almost identical frequencies, and the temperature at 
a fixed frequency changes by less than a factor of 1.5 
Therefore, the changes to these microphysical parametrs 
do not affect our earlier conclusions. From an astrophys- 
ical point of view, the most sensitive unknowns in this 
complicated calculation are the hyperon population (thus 
the equation of state and stellar mass) and superfluid 
bandgaps. At a very simple level this makes sense be¬ 
cause the hyperons appear only after a threshold density 
and the bandgaps appear in exponential factors rather 
than power laws. 
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IV. DISCUSSION 

We have extended previous investigations of r-modes 
in accreting hyperon stars in LMXBs as persistent 
sources of gravitational waves, focusing on the bulk vis¬ 
cosity which is needed to prevent thermal runaway. We 
have used improved microphysics compared to previous 
treatments, and have accounted for the most important 
macroscopic correction due to rotation of the star. 

We find that persistent gravitational wave emission is 
quite robust. Even the stiffest hyperonic equations of 
state in relativistic mean field theory produce enough 
damping to stop the runaway and persistently radiate, 
although some of them require neutron stars somewhat 
more massive than 1.4 Mq. Stars below about 1.3 Mq 
are not likely to be persistent sources regardless of the 
equation of state. The mass thresholds are somewhat 
more favorable for lower superfluid bandgaps than for 
higher bandgaps. Other details of the microphysics are 
found to be considerably less important. Our results 
seem robust for typical values of the crust-core viscosity 
and of saturation amplitude of r-modes, though this is¬ 


sue requires further investigation. Stars with high masses 
and stiff equations of state could exist in thermal and 
torque equilibrium at temperatures down to 10^ K. 

One possible avenue for substantial improvement is the 
hydrodynamics. We used a fluid expansion in the dis¬ 
sipation integrals which is a reasonable approximation 
for a Newtonian normal fluid, but not for a superfluid 
which can be much more complicated due to multiple 
components and entrainment. Work is underway to deal 
with the superfluid problem in general ([E3 and refer¬ 
ences therein). 
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